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A NOTE ON LAMBERT'S THEOREM 



This note combines all the various cases of Lambert’s Theorem into a 
single form which is particularly convenient for numerical work. This is made 
possible by appropriate choice of parameter and independent variable. 

Suppose a particle in a gravitational central force field has distances r t 
and r 2 at times 1 t and 1 2 from the center of attraction. Let c be the distance 
and 0 the central angle between the positions of the particle at the two times. 
Define 

s = (r t + r 2 + c)/2 
K = 1 - c/s 

q = ± K* 


The sign of q is taken care of by the angle 6 if we make use of 
e 2 = x | 2 + r 2 2 - 2r 1 r 2 cos 6 

to derive 

q = j^( r i r 2 ) W / s J cos ( 0 / 2 )- 

We further define 

G = universal gravitational constant 
M = mass of attracting body 
M = GM 

a = semimajor axis of transfer orbit 
E = -s/2 a for elliptic transfer 
= s/2 a for hyperbolic transfer 
T = (Syx/s)* (t 2 - tl )/s 
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m = number of complete circuits during transfer time. 


Note that - 1 < q < 1,0 < E <c0 for hyperbolic transfer, -1 _< E 
elliptic transfer, and E = 0 for parabolic transfer. Also 0 £ 0 £ v i 
and rr < 0 £ 2 77 if - 1 £ q < 0. 

Lambert’s Theorem 1 for elliptic transfer gives 
T = (~E )" 3/2 [ 2 m 77 + a - /3 ~ (sin a - sin/3)] 

E = ~ sin 2 (a/2) , 0 < a < 2tt 
sin(/3/2) = q sin(a/ 2 ), _ 77 £ /3 < 77 


For hyperbolic transfer, 

T - (E )' 3/2 [y - 8 - (sinh y ~ sinh 8)] 

E = sinh 2 (y/2) 
sinh ( 8/2 ) - q E 1/ 2 


If E is chosen 2 as the independent variable, a is ambiguous. We 
ambiguity by choosing as the independent variable 

x = cos (a/ 2 ), - l£x < l, 

- cosh ( 7 / 2 ) , x > 1 . 

For both elliptic and hyperbolic transfer 
E = x 2 - 1. 


< 0 f 71* 

0 1 q 1 1 


( 1 ) 


( 2 ) 


avoid any 
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For the elliptic case let 

y = sin (a/2) = (~E) H 

z = cos t/3/2) = (1 + KE)* 
v f = sin (54) (a- (3) = y(z - qx) 

g - cos (14) (a - /3) = xz - qE 

0 < a ~ < 2tt since 0 < f < 1 

h = (14) (sin a - sin /3) = y(x - qz) 

A. = arctan (f/g), 0 < k < v 
It then follows for the elliptic case that 

T = 2 (mn + k - h)/y 3 <3) 

For the hyperbolic case let 

y = sinh (y/2) = E* 4 

z = cosh (8/2) = (1 + KE)* 

f = sinh (14) (7-8) = y ( z - qx) 

g = cosh (14) (7-8) = xz - qE 

0 £ 7 - 8 < 00 since 0 £ f < 00 

h = (14) (sinh 7 - sinh 8) = y(x - qz) 

(14) (7-8) = arctanh (f/g) ' 

= (14) In [( f + g)/(g ~ f)] 

= (14) In jjf + e) 2 /(g 2 ~ f 2 )] 

= In (f + g) 
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Thus for the hyperbolic case 


T - 2 [h - In ( f + g)]/y 3 (4) 

When m - 0, equations (1), (2), (3) and (4) break down for x = 1 and suffer 
from a critical loss of significant digits in the neighborhood of x = 1 . To 
remedy this (1) is written 


T = 0(-E) ~ qK0(-KE), 

0(u) = 2jarcsin u'^ - u' A (l - u) 1 ^ . 

Replacing arcsin u A and (1 - u) M by series 3 , 


<fi( u) - 4/3 + ^ a„ u n , j u | < |, 

n= 1 


(5) 


a n = 1 • 3 • 5 ••• (2n - l) / /2 n_2 (2n + 3) n! 

A similar procedure produces the same series for the hyperbolic case. In 
fact (5) holds for the elliptic (m = 0) , parabolic, and hyperbolic cases provided 
0 < x < 2. 

It is now apparent that, given q and x , the following steps produce T for 
all cases: 


1. K = q 2 

2. E = x 2 - 1 

3. p = | E | 

4. If p is near 0 , compute T from (5). 

5. y - p' A 

6. z = (1 + KE)* 

7. f = y(z - qx) 
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8. g = xz - qE 

9. If E < 0, A. = arctan (f/g), d = mm + A., 0 < A. £ tt 

If E> 0, d = In ( f + g) 

10. T = 2(x - qz - d/y)/E 

The following formula for the derivative holds for all cases except for 
x = 0 with K = 1 and for x = 1. 

dT/dx = (4 - 4qKx/z - 3xT)/E 

If x is near 1 , the series representation should be differentiated. If q - 1 
we have a left-hand derivative of -8 and a right-hand derivative of 0 at x = 0 . 
If q = -1 we have a left-hand derivative of 0 and a right-hand derivative of -8 
at x = o. (See Figure 1. ) 

In the derivation of Lambert's Theorem for the elliptic case a and are 
defined in such a way that 

E 2 - Ej - a - y3 + 2m7T (5) 

where E x and E 2 are the values of the eccentric anomaly at times t x and t 2 . 
Thus from equation (1) 

E 2 - E x = (-E) 3/2 T + sin a - sin y3 

= y 3 T + 2y(x - qz). (6) 

We now obtain a formula for the scalar product 

Sj = r, • v, = rjVj sin «/» l 

of the position and velocity vectors at time t x , v x and 'p l being the speed and 
flight path angle. 
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Kepler's equation can be written 4 


(f^T (t 2 - t,) = E 2 - Ej * S,[l - cos(E 2 - E ,)]/(„.)' 

-(l - rj/a) sin (E 2 - Ej. 


Substituting a = - s/2E, t 2 - t t - s 3/2 T/(8 /x) % , and making use of (5) and 
(6) we have, after some algebra, 


St 


- (2 /j,s)' A qz(s - r x ) - x (s - 



A similar procedure produces the same formula for S t in the hyperbolic 
case. It also holds for the parabolic case. 


Figures 1 and 2 show T as a function of x for elliptic and hyperbolic 
transfer, the parabolic case occurring for x = 1. We suggest the reader com- 
pare these curves with those in Reference 2 showing T as a double-valued 
function of E with infinite slope at E = -1. 

No solutions of Lambert's equation exist in the shaded regions of figures 1 
and 2. x = l (m > 0) and x = -1 are vertical asymptotes. T 0 as x - 
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